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1. Introduction 

The QCD phase diagram and the location of a possible critical point (CP) are subject to intense 
theoretical and experimental research. Lattice QCD calculations support a crossover at vanishing 
baryon densities [1], while model studies indicate the existence of a first-order phase transition 
at high baryon densities ending at a CP [2]. In [3] it was shown how to detect the CP in heavy 
ion collisions by searching for divergences in event-by-event fluctuations of transverse momentum 
or particle multiplicity, an ansatz that has recently been refined to higher order cumulants [4, 5]. 
However, it is expected that finite size and time effects and possibly initial state fluctuations will 
crucially weaken the expected signals [6, 7]. On the other hand, the nonequilibrium evolution 
during a heavy-ion collision will enhance effects at the first-order phase transition, where spinodal 
instabilities may produce domain formation and clustering in energy and baryon density [8-10]. 
Hadronization of these clusters will lead to large non-statistical fluctuations in the hadron rapidity 
density within single events, providing an important observable signal for upcoming experiments 
at FAIR and NICA[11]. 

A successful dynamical model to study effects at the QCD phase transition in nonequilibrium 
has to go beyond usual hydrodynamics that includes the phase transition only in the equation of 
state [12]. A novel approach that includes the dynamics of the order parameters explicitly is given 
by chiral fluid dynamics [13-17]. We recently extended this model with the Polyakov loop to 
consider both the chiral and the deconfinement transition [10, 18]. 

2. Polyakov-chiral fluid dynamics (P;^ FD) 

The basic idea of the model is to explicitly propagate the sigma field and an effective Polyakov 
loop as the order parameters of the chiral and deconfinement phase transition. A fluid dynamically 
expanding medium of quarks and antiquarks provides the locally thermalized background for these 
fields. This enables us to study relevant effects at the CP and first-order transition in a dynamical 
system of finite size. 

We use the Polyakov loop extended quark meson model [19] with the Lagrangian 

^^q[i (r^a^ - ig.fAo) -ga]q+^{d^Gy-U (a) - ^ {£, I) , (2. 1) 

where q = {u,d) is the constituent quark field, Aq the temporal component of the color gauge field, 
O" the mesonic field and i the Polyakov loop. The pion degrees of freedom are neglected throughout 
this work. The potential for the sigma field is the usual "Mexican hat" 



U{a) = 




(2.2) 



and the temperature dependent Polyakov loop potential is chosen in a polynomial form [19, 20]: 




(2.3) 
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Integrating out the quark degrees of freedom in the partition function ^ gives us the effective 
potential: 



V 



ln^ = C/ + ^ + ^ 



qq 



(2.4) 



Here, the quark contribution Qqq determines the local equilibrium pressure of the quark fluid. In 
mean-field approximation and at zero chemical potential it reads [19]: 



= -4NfT 



J {In 



In 



l + 3^e-^^ + 3^e-^^^ + e-^^^ 



(2.5) 



We tune the strength of the transition by varying the quark-meson coupling g. This allows 
us to study first-order phase transitions and transitions through the CP at vanishing baryochemical 
potential. Figure 1 shows the effective potential for g = 4.7 (first-order) and g = 3.52 (CP) at the 
respective transition temperature (cf. Ref. [10]). Note that in general one has to choose g such that 
the product go resembles the constituent quark mass in vacuum, leading to a value of g ^ 3.3. 




Figure 1: la Effective potential for g = 4.7, corresponding to a first-order phase transition at Tc = 
172.9 MeV. lb Effective potential for g = 3.52, corresponding to a CP scenario at Tc = 180.5 MeV. Both 
figures are adopted from [10]. 



One can quantify this behavior by calculating the chiral and Polyakov loop susceptibilities Xao 
and Xii' In Fig- 2 they are shown for three different couplings. We find divergent susceptibilities 
for g = 3.52 indicating a chiral and deconfinement CP. 

Within the two-particle irreducible effective action formalism we self-consistently derived the 
coupled dynamics for the sigma field and the quark heat bath [15]. We obtained a Langevin equa- 
tion for the sigma field with temperature dependent damping rja and stochastic noise term £,a that 
are connected via a dissipation fluctuation relation 

d^d^'a + ria{mG+^ = , (2.6) 

OG 

{Ut,x)Ut',^)) = ^5{t-t')5ix-^)mariaC0th{^) . (2.7) 
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Figure 2: 2a Chiral susceptibility as a function of temperature for different coupling strengths. 2b Polyakov 
loop susceptibility as a function of temperature for different coupling strengths. 



For the Polyakov loop field we deploy a relaxation equation which also contains stochastic noise: 

^,ddT' + ^^^,^ (2.8) 

{^e{t,xMt',^))T^ = Uit-t')5ix-^)2l^eT . (2.9) 

Note here that the Polyakov loop is originally defined only in equilibrium and it is not a priori clear 
what the correct dynamics are [20]. This approach is therefore purely phenomenological. A similar 
ansatz with an additional kinetic term has been pursued in [21, 22]. The damping coefficient rji 
is set to a value of 5 / fm. Results are sensitive to this choice only in the vicinity of the first-order 
transition temperature [10]. 

The quarks are propagated via the equations of ideal relativistic fluid dynamics: 

d^T^^'' = Sl + Sj , (2.10) 

with source terms 5^ and SJ describing the energy transfer from the fields to the fluid via damping. 
The energy transfer due to stochastic fluctuations needs to be estimated numerically [10, 17]. 



3. Numerical results 



3.1 Equilibration in a box 

We study several temperature quenches in a cubic box with periodic boundary conditions. 
Both fields are initialized at some global Ti^i > Tc, with Tc being the respective critical temperature. 
Then the temperature is quenched to a value T < Tc and the energy density and pressure of the 
quark fluid are calculated. We let the coupled system evolve and relax. As pressure gradients are 
small within this setup, we expect the dynamics to be dominated by the fields. The solid red curves 
in Fig. 3 show the volume and event averaged sigma fiel (a) for equilibration near the transition 
point for both first-order and CP scenarios. At the first-order transition the significant delay in the 
relaxation time is caused by the large barrier separating the degenerate minima. Critical slowing 
down can be observed near the CP, where the vanishing of T]cr causes oscillations and prevents the 
field from relaxing to its equilibrium state. Similar effects occur in the Polyakov loop [10]. 
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Figure 3: 3a Equilibration of the sigma field for several quench temperatures T <Tc through the first-order 
transition. The barrier between the minima in the potential increases the relaxation time when the system 
relaxes near Tc = 172.9 MeV. We chose Tini = 180 MeV. 3b Equilibration of the sigma field for several 
quench temperatures T < Tc through the CP. Critical slowing down delays the dynamics and causes oscil- 
lations around the flat minimum when the system relaxes near Tc = 180.5 MeV. We chose Tini = 186 MeV. 
Both figures are adopted from [10]. 



Another critical phenomenon can be observed by studying the intensity of field fluctuations. 
These are given for the sigma and Polyakov loop field as [10, 23]: 

dNg ^ co^\5ak\^^\dtGk\^ dN^ ^ + 

d^k {2ny2cOk ' d^k {27if2(Ok ' ^ ' ^ 

Here 8ok and dtO^ are the ^h Fourier modes of 5(7 = <7 — Ceq and dfO and (O^ is the corresponding 
energy. We compare intensity histograms in the late stage of the evolution in the CP and first-order 
scenario in Fig. 4. For both order parameter fields we find a strong enhancement of long- wavelength 
modes at the CP compared to an equilibration near the first-order transition point. 




I A: I (GeV) |A:| (GeV) 

(a) (b) 

Figure 4: 4a Intensity of sigma fluctuations after equilibration att = 24 fm. In the CP scenario we find an 
enhancement of the soft modes. 4b Intensity of Polyakov loop fluctuations after equilibration at ^ = 24 fm. 
In the CP scenario we find an enhancement of the soft modes. Both figures are adopted from [10]. 



3.2 Fluid dynamic expansion 

To explore the influence of the expansion on the dynamics of the fields, an ellipsoidal region 
with a temperature T = 200 IMeV, above both transition temperatures, is provided as initial state 
of a fluid dynamic expansion. This is to resemble the situation after the collision of two heavy 
nuclei. Fields and fluid are again set to their respective equilibrium values and the system evolves 
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according to full (3+l)-dimensional fluid dynamics. During the expansion we observe supercool- 
ing and reheating in the first-order transition scenario. This supercooling causes an enhancement 

of nonequilibrium fluctuations (Ao") = ^J {[a — cJeq)^) and (A^) = ^J {{l — ltqf) in both order pa- 
rameters at the first-order phase transition, see Fig. 5. The second bump in the fluctuation strength 
near = 6 fm arises when parts of the system cross the transition temperature a second time after 
reheating. 
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Figure 5: Nonequilibrium fluctuations of the sigma field 5 a and Polyakov loop 5b are enhanced at the 
first-order transition compared to the CP scenario. 



3.3 Domain formation at the first-order phase transition 

We now focus on the evolution of a single event to learn more about the transition processes. 
To achieve this we introduce spatial correlations for the stochastic noise fields over volumes 
and to obtain a more physical behavior of these fluctuations. 

We show a slice in the transversal z = plane for the sigma field, Polyakov loop and energy 
density in Figs. 6, 7 and 8, each for early, intermediate and late times in the evolution. We see 
in the order parameters domains of the high- and low-temperature phases coexisting during the 
transition process. This phenomenon is typical for the first-order phase transition and does not 
occur in evolutions through the CP. It can be best observed in the sigma field, but also the Polyakov 
loop exhibits a bumpy structure during its evolution. This structure then translates to the energy 
density, leading to a significant amount of inhomogeneity and clumping. 
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Figure 7: Polyakov loop field in the z = plane for r = 1 fm 7a, ^ = 4 fm 7b, and ^ = 7 fm 7c during a 
first-order phase transition. Fig. 7b adopted from [10]. 
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Figure 8: Energy density in the z = plane for ^ = 1 fm 8a, ^ = 4 fm 8b, and r = 7 fm 8c during a first-order 
phase transition. Fig. 8b adopted from [10]. 



We expect this effect to become even stronger when we go to systems at finite baryon density. 
This would then provide an important experimental signal for the QCD phase transition, e. g. in 
non-monotonic multiplicity fluctuations of hadrons. 

4. Conclusions 

We presented the extension of nonequilibrium chiral fluid dynamics with a Polyakov loop 
to include effects of the deconfinement phase transition of QCD. We were able to observe typi- 
cal critical phenomena like critical slowing down and the enhancement of soft modes for systems 
equilibrating near the CP. For an expanding system cooling through the first-order phase transition 
we found evidence for the formation of a supercooled phase leading to subsequent reheating of the 
fluid. As a result, large nonequilibrium fluctuations evolve. For single events, we find significant 
difference in the evolution of fields and fluid between the CP and the first-order scenario. The latter 
one proceeds through the formation of domains in the order parameter fields leading to irregular- 
ities in the energy density. As a next step we investigate this effect for systems at finite chemical 
potential to provide relevant signals of the QCD phase transition for upcoming experiments at 
FAIR. 
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